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Abstract 

Buckling  is  a  significant  cause  of  structural  failure.   A 
theoretical  investigation  has  been  undertaken  into  the  dynamic 
instability  of  complete  spherical  shells  which  are  loaded 
impulsively  and  made  from  a  viscoplastic  material. 

The  critical  mode  numbers  obtained  in  this  thesis  are 
similar  to  those  obtained  for  cylindrical  shells  having  the  same 
R/h  ratios  and  material  parameters  while  the  threshold  velocities 
were  for  higher  viscoplastic  spherical  shells  than  for  rigid- 
plastic  spherical  shells. 
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Introduction 

The  analysis  of  dynamic  buckling,  of  elastic-plastic 
and  rigid  plastic  cylinders  has  been  investigated  by 
Abrahamson  and  Goodier  [1]  extended  by  Goodier  [2] ,  Goodier 
and  Florence  [3]  and  others  [4-6].   Florence  [7-8]  and 
Wojewodzki  [9-11]  extended  the  analysis  to  viscoplastic 
cylindrical  shells.   V7ork  on  spherical  shells  has  progressed 
through  elastic-plastic  and  rigid  plastic  materials  [12-14]. 
The  analytical  procedures  of  the  above  references  were  derived 
by  assuming  that  buckling  stemmed  from  the  growth  of  small 
imperfections  in  otherwise  uniform  initial  displacements 
and  velocity  fields. 

The  dynamic  buckling  of  a  thin-walled  complete  spherical 
shell  which  is  made  from  a  viscoplastic  material  is 
investigated  herein.   A  theoretical  procedure  is  developed 
by  first  considering  the  unperturbed  motion  arising  from 
a  spherically  symmetric  initial  velocity  field.   An 
examination  is  then  made  of  the  influence  of  small  perturbations 
in  this  otherwise  uniform  velocity  field  and  of  the  effect 
of  small  deivations  from  the  initial  sphericity  of  the  shell. 
It  transpires  that  certain  harmonics  grow  very  quickly  and 
cause  a  spherical  shell  to  exhibit  a  characteristic  wrinkled 
shape.   It  is  assumed  that  the  strains  which  are  associated 
with  the  perturbed  state  are  much  smaller  than  those 
associated  with  the  predominant  motion.   Thus,  no  unloading 
of  the  shell  material  occurs  before  the  time  when  all  the 
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initial  kinetic  energy  is  absorbed  by  the  spherical  shell  in 
the  predominant  mode.   Moreover,  it  is  further  assumed  that 
the  wrinkled  pattern  in  the  shells  is  established  prior  to 
this  time  and  is  not  modified  by  subsequent  unloading. 
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Chapter  2 
Basic  Equations 

The  equilibrium  equations,  compatibility  relations  and 
strain  and  curvature  relations  have  been  presented  in  refs.  [13] 
and  [14]  for  a  spherical  shell.   These  results  are  repeated  in 
this  section  for  convenience. 

It  may  be  shown  for  the  particular  case  of  a  thin-walled 
spherical  shell  that  the  strain  and  curvature  relations  which 
were  developed  in  reference  [12]  for  an  arbitrarily  shaped 
shell,  can  be  written 

%  =  I  (%  "  "  +  !%%>  (la) 

eg  =  =■  (cosec<}>  V,  Q  +  cot<f>  U  -  W  +  —  cosec2<J>  W,gW,Q)    (lb) 

e.Q  =  ^[cosectf)  u,Q  -  cotcj>  V  +  V,  ,  + 

|  coseccf.  (Wf(j)W/0  +  WfQWf  )]  (lc) 

%  -  W   W  +  %♦>  (ld) 

•  1  •  •  • 

<fl    =  —t    (W   +   cosec2*   W/flB   +   cot<|>  W,     )  (le) 


R 


and      k.q    =  —g-    [coseccj)  W,  ,„    -    cottj)   cosecij)   W'ol  (If) 

P.  ' 

when  the  displacements  are  defined  as  shown  in  Figure  1.   A 
straightforward  modification  of  the  general  equations  in 
ref.  [12]  for  a  shell  of  thickness  h  and  density  p  gives  the 
equilibrium  equations 
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(b) 


Figure  1:   (a)  Displacement  and  (b)  Membrane  Forces,  Bendinq 
Moments  and  External  Loads  on  a  Spherical  Shell. 
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V*  +  C°t4>  (N*  "  N9)  +  COSec<})  N*6'0  " 

phRU  +  P1R  =  0  (2a) 


cosec*  Ne,e  +  N^,^  +  2  cot*  N^g  - 

phRV  +  P2R  =  0  (2b) 


and        R(N^  +  NQ)  +  N^W,^  +  cot*  W/(j))  + 

Nq  cosec2*  W,qq  + 
cosec*  N^g  (Wf(J)g  +  W,e(J))  + 

(V*  %  +  W'6N(f.e'<(.  cosec*  + 

w'*  N*e'e  cosec*  + 

cosec2*  Ne,e  W,e)  -  (M^^  + 

2  cot*  M^  -  cot*  Mg,^  + 

2Mg  +  cosec2*  Mg,gQ  + 

2  cosec*  M^g^g  + 

2  cosec*  cot*  M.0/q)  -  R2 phw  + 

R2P3  =  0  (2c) 

which  are  consistent  with  equations  (la)  and  (If)  according  to 
the  principle  of  virtual  work  provided  the  influence  of 
transverse  shear  deformation  is  disregarded. 

The  strain  and  change  of  curvature  relations  (la)  to  (If) 
satisfy  the  compatibility  equation 
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ee'<Jxj)  +  cosec2*  Vee  +  2  cot*  e9'4,  " 
cotD  e^  +  2e^  -  coseccf)  e^,^  - 

•  •  • 

coseccj)  cotij)  e,„,„  +  R(<A  +  <fl)  =  0  (3) 

for  infinitesimal  displacements  [12] . 

The  constitutive  equations  for  spherical  shells  which  are 

made  from  rigid  viscoplastic  material  are  developed  in  Appendix  A. 

Equations  (1.8a)  to  (1.3f)  for  a  spherical  shell  made  from  rigid 

viscoplastic  material  can  be  written  in  the  form 

kh 

N<j>  "  —  (2e<j>  +  ee)  (4a) 

kh 

Ne  =  T  (2ee  +  V  <4b> 

Y 

<h  eAO 

(4c) 


4)6         Y*       2 

kh3 

j 

M        —                         f°i' 

4-      k"      1 

*         12Y*            * 

+      <eJ 

kh3 

. 

M          —                                f  T|/- 

+     ^      ^ 

6             TO      *                 9 

12y* 

+     K<t>' 

and 

kh'      K<$>6 
M        — 

*9         12Y*       2 

where 

k  =   a   //3 

and 

Y*    =   Y(l    -    k//j7)  . 

(4d) 


(4e) 


(4f) 

(5a) 

(5b) 

The  method  of  solution  used  here  is  similar  to  that  used  in 
refs.  [12]  and  [13].   It  is  assumed  that  the  displacement  in  a 
complete  spherical  shell  consists  of  dominant  parts 
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(W,  U  =  V  =  0)  and  perturbation  parts  (w' ,  U' ,  V)  so  that 

W(4»f9fT)  =  W(T)  +  W  (<|>,efT)  (6a) 

U(4>,e,x)  =  V  (4>,0,t)  (6b) 

and        V(<J>,8,t)  =  V  (*,  6 ,  t)  .  (6c) 

Thus,  the  dominant  strains  and  curvature  changes  are 
_    w      _ 

e4)  =  ee  =  "  r  '   e<(,e  =  °'  (7a_c) 

and       k^   =  kq   =  iJr  ,     ^0  =  0  (7d-f) 

respectively,  while  the  corresponding  perturbation  quantities 
are  given  by  equation  (la)  to  (If)  with  primes.   Similarly, 

Ni(<i>,e,T)  =  N(t)  +  Nj(*f6fT)  (8a) 

Ne(<f),e,T)  =  N(T)  +  N£(<J>,9,t)  (8b) 

N({)e(4),e,T)  =  »je(*,efT)  (8c) 

M^(*,8,T)  =  M(T)  +  Ml($,6fT)  (8d) 

Me(4>,e,T)  =  M(T)  +  M'(4>,9,x)  (8e) 

and       M^e(#f6fT)  =  M^Q  (<J>,  6,t)  .  (8f) 

If  equations  (6)  and  (8)  are  substituted  into  equation  (2c), 
then  the  dominant  and  perturbation  transverse  equilibrium 
equations  are 

N      M  „ 

2 2  —  +  P-  -  phW  =  0  (9a) 

R      R2     3 


-  16  - 


and  R(N'    +   N^)    +   N(W'     .    +   cot<|>  W)  .    .+   cosec2<J>   W"     Q)    - 

tMi; <h  +  2  cot*  Mi'd,  ~  cot*  Me'0  + 

2M»  +  cosec2cJ>  M'QtQQ   +  2  cosec*  M^,^  + 

2  cosecij)  cot(J)  M'  fQ]  -  phR2W  =  0  (9b) 

respectively,  when  neglecting  higher  order  quantities. 

It  may  be  shown  that  in-plane  equilibrium  equations  (2a) 
and  (2b)  are  satisfied  identically  when 

N*   =   R2"    [F   +    oose°2^   F'<j><t>   +   cot(i>   F'<j)]  (10a) 

Ne  =  ^   [F  +  F'^]  (10b) 

and       N.q  = J"  coseccj)  [F,  ,  Q  -  cot<J>  F,Q]  (10c) 

where  F  is  a  stress  function  and  U=V=P,  =P2=0.   Now 
substituting  the  inversions  of  equations  (4a)  to  (4c)  and 
equations  (12a)  to  (12c)  and  (Id)  to  (le)  into  the  compatibility 
equation  (3)  gives 

3(V2  +  2)u,T  +  [2(V2  +  1)  -  1] (V2  +  2)F  =  0        (11a) 

where      u  =  —     F  =  -^ C2  =  —     t  =  —  t  (llb-e) 

wnere      u    R  ,  i         khRC  ,   ^     p  '       R 

and       V2  (  )  =  (  )     +  cot<M  )  ,   +  cosec2  ({.(),  Qe .         (llf) 

The  transverse  equilibrium  equations  (9a)  and  (9b)  can  be 
written  in  the  form 

~uo'tt  +  2(1  +  a2)°  +  P  =  0  (12a) 
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and 


u,    -  B2 (V2  +  2)F  -  aV2u  + 

XT 


cT3  [2(V  +  1)  +  1]  (V  +  2)u,   =  0 


(12b) 


respectively,  when  using  equations  (1),  (4),  (6)- (8)  and  (10) 
and  where 


and 


u   = 
o 


32  = 


w 

R    ' 

a    -     N 

kh    ' 

PoR 
P   =   -=^_    , 

kh 

2                1 

a      =  — 
12 

R 

>.         J 

C 

* 

Ry' 


(12c-g) 
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Chapter  3 
Static  Behavior 


Equation  (12a)  predicts 


a   =   -   — 
2 


(13) 


when  a2  <<  1,  while  equations  (11a)  and  (12b)  become 


[2(X   -  1)  +  1]  (A   -  2)c 
1  v  n         J   n      mn 


3(A   -  2)a   ,   =  0 

n      ran '  t 


(14a) 


and 


32  (A   -  2)c    +  5A   a    + 
n      mn     n  mn 


a2B2  [2(A   -  1)  -1]  (A   -  2)a   ,   =  0 
n  n      mn ' t 


(14b) 


respectively,  if  it  is  assumed  that 


u  =  aon  Vcos<}))  + 


jn 


P  (cos<J>)  (a    cos  m6  +  b    sin  m6) 
n        mn  mn 


(15a) 


and 


F  =  c   P  (coscf>)  + 
on  n     r 


P  (cos*) (c    cos  m9  +  d    sin  m6) 
n     Y    mn  mn 


(15b) 


where 


A   =  n(n  +  1) 
n 


(15c) 


and  P   and  P   are  Legendre  and  associated  Legendre  polynomials 

of  degree  n  and  order  m.   Two  additional  equations  are  also 

obtained  but  these  are  identical  to  equations  (14a)  and  (14b) 

provide  b    and  d    are  substituted  for  a  „  and  c   .   If  c  „ 
c  mn      mn  mn      mn       mn 

is  eliminated  from  equations  (14a)  and  (14b) ,  then 
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3  to   a   ,   +  oX      a 
n  mn '  t     n  mn 


(16a) 


where 


3(X   -  2) 

to   = +  a2(2A   -  3)  (A 

n     2X      _ 

n 


2). 


(16b) 


Now,  substituting  equation  (13)  into  (16a)  and  integrating 
gives 


fP  xn    1  C  .  1 

2  ,7,_  B2  R 


w 


a  .  =  e 
mn 


n 


+  K 


(17) 


where  K  is  a  constant  of  integration  which  is  related  to  the 

7  C 

initial  imperfections  in  a  shell.   3   which  equals    %  depends 
on  strain  rate  as  y*  depends  on  strain  rate  as  shown  in 
appendix  B.   With  a  strain  rate  near  zero,  as  would  be  the 
case  for  static  buckling  32  goes  to  infinity  but  1/3 2  becomes 


Ry*    R 
C     C 


e   y 

e  ' 


e   +  y 
e    ' 


(18) 


In  the  limit  of  e      very  small,  e   <<  y 


limit 

e  ->  0 
e 


(19) 


where    e      is   the    initial    strain   rate, 
e 


becomes 


a™*,    =   e 

mn 


|  ^  e   tl 

2    wn      ej+K 


Therefore  equation  (17) 


(20) 


when  e   <<  y 
e 


This  equation  is  of  the  form  of  creep  buckling. 
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Chapter  4 
Dynamic  Buckling 


4.1   Dominant  Equation 

From  equation  (12a) 


u 


o'  TT 


2a 


(21) 


when  a2  <<  1  and  P  =  0.   a  is  not  time  invariant  and  the 
solution  of  the  above  equation  will  be  done  in  two  parts. 
The  first  part  will  assume  the  value  of  a  to  be  constant  with 
time  and  the  second  will  give  the  exact  solution,  then  the  two 
solutions  will  be  equated  at  the  final  displacements.   For 
buckling  the  dominant  stress  is  defined  as 


N  =  -  ha. 


and  equation  (21)  becomes 


(22) 


U  ,    =  -  2/3 

O  TT 


(23) 


The  solution  of  equation  (23)  is 


V 


u   = 
o 


1  - 


2t, 


(24a) 


when  satisfying  the  initial  conditions  W 
t  =  0  and  where 


0  and  W  =  V   at 
o 


V 


2/3C 


r         \ 

a 
o 

a 

(24b) 


is  the  time  when  W  =  0, 
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The  exact  solution  uses  equation  (4a) ,  (4b) ,  (7a)  and  (7b) 
to  revise  a   to  reflect  the  time  dependence  as  follows 


N 


kh 

—   3 
* 

Y 


W 
R 


(25) 


Solving  for  a  and  using  equations  (12c)  and  (2.9)  yields 


a   = 


N 
kh 


3    _  _C_ 

2    J  Ry  Uo't 


.where 


2  v:  u  ,  . 
R   o  T 


(26) 
(27) 


Substituting  equation  (26)   into  (21)  produces 

6C 

u  ,    +  —  u  ,   =  -  3 

O  TT     „     O  T 

Ry 


(28) 


the  equation  for  the  exact  solution.   The  solution  of  equation 
(28)  is 


o    12 


ry 


2V 


1 


1  + 


RY  J 


1  -  e 


6C 

"r7t 


ry 


2C 


(29a) 


which  satisfies  the  initial  conditions  W  =  0  and  W  =  V  at  t  =  0 

o 


and 


RY 


Tf    6C 


log. 


2V 


l  + 


RY  J 


(29b) 


is  the  time  when  W  =  0.   Setting  u  (tf)  for  the  first  solution 
equal  to  u  (^f)  for  the  exact  solution  and  solving  for  o   /o 
one  obtains 


/3 


RY 


V 


2V 


ry 


-  log 


1  + 


2V  1 
o 

ry 


(30) 


-  22  - 


therefore,  the  dominant  solution  may  be  taken  as  equation  (24a) 
with  (24b)  and  (30)  used  with  it. 
4.2   Perturbation  Equation 

Using  equations  (15a-c)  equation  (12b)  becomes 

a   ,    -  B2  (A   -  2)c    +  —  A   a    + 
mn  tt       n      mn    ,  ,   n   mn 

kh 

a2B2[2(X   -  1)  -  1]  (A   -  2)a   ,   =  0    (31) 
n         J   n      mn'  x     .' 

Solving  equation  (14a)  for  c    and  substituting  into 
equation  (31)  produces 


a   ,    +  3  (jo   a   ,   +  —  A   a    =    0 
mn'xT      n   mn'x    .  .   n   mn 

kh 


(32a) 


where 


X      =   n(n  +  1) 


n 


(32b) 


3(A   -  2) 


u 


n 


n 


2A   -  1 


+  a2  (2A   - 


n 


3)  (A   -  2) 
n 


n 


(32c) 


RY" 


1      C 

+  — 


2U0'T     ^ 


(32d) 


With  a  change  of  variable  equation  (32a)  becomes 


r 


"mn'?? 


T_  0) 

f   n 


C      C 

+  — 


2V0S    RY. 


1     i   r 

mn'  K 


3    3V 

- .+  — 2 

2    Ry 


A   a 
n  mn 


(33a) 


where 


5  =  1  -  t/t 


(33b) 
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N 
kh 


—  + 
2 


3V 


RY 


(33c) 


and 


V  ■ 


(33d) 


Putting  equation  (33a)  into  a  format  similar  to  Florence  [8] , 
one  gets 


amn'CC 


O 


(2) 


o(1)  +'2L 

~n        r 


a   ,  c    -    R   a    =0 
mn  t,  n  mn 


(34a) 


where 


G(1)  =  T,  ^^ 

"n       f  Ry   n 


(34b) 


Q 


(2) 


n 


=  T 


C   - 

j.   to 

f  2V„   n 


(34c) 


and 


R2  =  TJ 

n     f 


3  ,  3VP 
2     Ry 


n 


(34d) 


As  R   is  still  a  function  of  time  (£) ,  in  order  to  solve  the 

(3    3vo   1 

equation  using  Florence's  solution,  the  value  of  -j  +     R   5   was 

set  as  a  constant  with  respect  to  £.   A  value  of  3/2  was  used  as 
we  are  primarily  interested  near  £  =  0  which  negates  the 
influence  of  the  second  term,  thus 


R2  =  3/2  t2  X^ 
n    '    f  n 


(34e) 


As  given  in  Florence  [8]  the  solution  to  equation  (34a)  is 


a(£)  =  A(£)  a   +  Bm  (?)Bn 
mn       mn     n    mn    n 


(35) 
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where     ^(O  - 


Til   +   a„  +  b  ) 

n     n 

r(2  +  b  ) 
n 


exp{-  (A   +  Q    )    +    A  O 
*  n    n     n 


1+b 


n 


n 


-  {_(An  +  1  +  bn)  U(fin)  + 


a   u*  (n  )}m(q  ?)  +  {(A  +  1  +  b  )  M(n  )  + 

n    n     n        n        n     n 


fi   M'  (ft  )  }  U(ft  £) 
n     n      n 


(36a) 


and 


B   (?) 
mn 


r(l  +  a   +  b  ) 

__ n_   n 

T(2  +  b  ) 
n 


exp{- (A   +  Qn)  + 


V> 


1+b 


(ft  K) 

n^ 


n 


U(ft  )M(ft  £) 
n     n 


M(fi  )u(ft  K) 

n  -    n 


2t 


(36b) 


In  equations  (36a)  and  (36b),  M(fi  ?)  and  U(ft  £)  are  shorthand 
notations  to  represent  Rummer  functions  [15-17]  (confuent 
hypergeometric  functions)  of  the  first  and  second  kind, 
respectively,  that  is 

M(ft  ?)  =  M(l  +  a   +  b  ,  2  +  b  ,  ft  £)        (37a) 
n^  n    n       n    n 

U(ft  £)  =  U  (1  +  a   +  b  ,  2  +  b  ,  ft  £)•       (37b) 
n  n    n'       n    n 

Similarly 


M(nn)  =  M(l  +  an  +  bn,  2  +  bn,  Qn) 


U(J2  )  =  U(l  +  a   +  b  ,  2  +  b  ,  ft  ) 
n  n    n'       n    n 


(37c) 
(37d) 
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and  primes  over  M(Q  )  and  U(ft  )  denote  differential  coefficients 

n         n 

with  respect  to  R  C  evaluated  at  5  =  1 .   The  remaining  symbols 


are 


Q 


(2) 


a   = 

n 


"n 


0 


(1) 


n 


-  1 


n 


n 


b  =  o<2> 

n    n 


(37e-f) 


n  = 
n 


Q(1)   +  4R2 
n        n 


1/2 


0 


(1) 


0. 


A  = 
n 


n 


n 


(37g-h) 
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Chapter  5 
Results  and  Discussion 

Unfortunately,  no  experimental  investigations  appear  to 
have  been  published  on  the  dynamic  instability  of  complete 
spherical  shells  and  until  reciently  no  theoretical  work  had 
been  published.   The  results  reported  herein  will  then  be 
compared  with  previous  experimental  investigations  into  the 
behavior  of  impulsively  loaded  cylindrical  shells  and  in  theory 
with  the  publications  now  in  press.   Equation  (5.2)  from 
Appendix  E  predicts  mode  numbers  which  are  3.173/3.624  =  0.876 
times  as  large  as  Jones  and  Ann's  [13]  estimates  for  similar 
aluminum  spherical  shells  and  are  3.173/74'  =  1.08  times  larger 
than  Florence  and  Vaughan's  [4]  estimates  for  similar  aluminum 
cylindrical  shells.   However,  as  shown  in  Figure  2,  a  plot  of 
the  first  five  sets  of  data  from  Table  1,  the  equation  for  the 
critical  harmonic  is 

n  =  2.0(R/h) *605  (38) 

which  compares  favorably  with  Lyons  [6]  curve  with  a  slope  of 
0.625. 


R/h         Critical  Harmonic  n 

10  8 

18.25  12 

30  15 

50  21 

100  32 

Table  1:   A  Comparison  of  the  Critical  Mode  Numbers  for  a  Visco- 
plastic  Spherical  Shell  with  Differing  R/h  values,   _, 
with  R  =  1.460  in.,  V0  =  9357.6  in/sec,  Y  =  6500  sec" 
a      =  46,000  psi  and  p  =  0.0002523  lb-sec2/in.  "•  . 
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Figure  2:   Critical  Mode  of  Buckling  as  a  Function 


of  R/h  with  VQ  =  9357.58  in/sec  and 


Y  =  6500  sec 


-1 
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Figure  3  shows  the  growth  of  the  amplitude  function  with  regard 
to  time  for  the  continuous  spectrum  of  the  mode  numbers  n. 
Showing  addition  figures  later;  this  growth"  is  faster  for 
larger  V  ,  larger  y   and  for  later  stages  of  the  process 
(small  K) •   The  change  of  preferred  buckling  mode  is  shown 
by  a  dashed  line.   It  is  evident  from  Table  2a  that  equation  (38) 
provides  a  reasonable  estimate  of  the  preferred  mode  of  the 
amplification  function  which  is  obtained  numerically  from 
the  corresponding  theoretical  analysis  using  Kummer ' s  functions 
and  numerical  integration  of  the  equations.   The  use  of  the 
Kummer ' s  functions  for  solution  ran  into  problems  because  of 
the  size  of  the  functions  and  that  differences  of  these  very 
large  numbers  were   removing  the  first  five  or  six  significant 
figures.   Even  double  precision  was  not  eliminating  the 
problem,  for  the  asymptotic  expansion  [16]  is  a  truncating 
series  and  even  with  double  precision  can't  give  the  necessary 
significant  figures. 

The  Florence  method  gave  good  results  that  correlated 
well,  up  to  four  significant  figures,  with  the  numerical 
method.   The  problem  with  the  Florence  method  occurred  when 
the  values  of  inputs  to  the  Kummer  function  became  so  large 
i.e.  greater  than  8  that  one  has  to  use  an  asymptotic 
expansion.   Also  the  last  couple  of  computations  with  values 
just  less  than  8  one  loses  some  accuracy.   As  was  observed 
by  Florence  [8]  and  Wojewodzki  [11]  the  critical  modes  of 
viscoplastic  cylinders  were  constantly  lower  than  modes  of 
rigid-plastic  cylinders.   Likewise  as  shown  from  Table  2a, 
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Material 


Ref . 


(ksi) 


(lb-sec2/in. * ) 


sec 


-1 


Steel  1015 

[1] 

102 

AL  2024  T3 

[1] 

56 

AL  6061  16 

[1] 

44 

AL  60  61  T6 

[5] 

44 

AL  6061  T6 

[6] 

46 

Steel  1015 

[7] 

30 

0007324 

40.4 

0002523 

6500 

0002523 

6500 

0002525 

6500 

0002523 

6500 

0007324 

40.4 

Table  2b:   Values  of  a  ,  p  and  y   that  Correspond  to  Materials 
in  Table  2a. 
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Material 
No. 


EO. (24b) 


C(103in/sec) 
EO. (lid) 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 


094 
110 
108 
085 
158 
140 
182 
146 
153 
166 
173 
171 
097 
105 
127 
130 
142 
153 
164 
105 


16.1 
16.3 
16.3 
16.0 
15.4 
15.3 
15.5 
13.5 
13.5 
13.6 
13.6 
13.6 
13.2 
13.2 
1-3.4 
13.4 
13.8 
13.8 
13.9 
13.6 


Table  2c:   Values  of  Tf  and  C  that  Correspond  to  Materials 


in  Table  2a.    (Continued) 
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Material  Tf 

No.  EO.'(24b) 


C(103in/sec) 

EO. (lid) 

13 

7 

13 

.7 

14 

.1 

14 

0 

13 

9 

14 

0 

13 

7 

13 

9 

14 

1 

8 

5 

8 

4 

8 

6 

8. 

3 

8. 

2 

8 

7 

8. 

5 

8. 

5 

8. 

7 

8. 

4 

8. 

3 

8. 

2 

21  .120 

22  .128 

23  .200 

24  .179 

25  .169 

26  .179 

27  .123 

28  .156 

29  .198 

30  .118 

31  .103 

32  .145 

33  .082 

34  .071 

35  .157 

36  .119 

37  .115 

38  .158 

39  .103 

40  .087 

41  .076 


Table  2c:   Values  of  f   and  C  that  Correspond  to  Materials 
in  Table  2a. 
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the  critical  modes  of  spherical  shells  are  considerably  lower 
than  those  observed  experimentally  on  cylindrical  shells  with 
the  same  material  parameters  and  R/h  ratios  which  were  tested 
by  the  authors  of  references  [1,  5,  6,  7],  with  the  further 
observation  that  correlation  seems  to  be  closer  at  lower  R/h 
values.   Of  course,  n,  in  the  present  case,  refers  to  the 
critical  or  most  amplified  harmonic   in  the  displacement 
field  of  a  complete  spherical  shell  which  is  expressed  in  terms 
of  Legendre  and  associated  Legendre  polynomials  (equation  (15a) ) , 
while  in  reference  [4]  on  cylindrical  shells,  it  is  associated 
with  trignometric  functions. 

The  concept  of  threshold  or  cirtical  impulse  has  been 
used  by  numerous  investigators  to  mark  the  smallest  impulse 
that  a  structure  can  tolerate  without  excessive  deformation. 
It  is  evident  from  Figure  4  which  is  a  plot  of  Table  3  that 
the  displacement  amplification  function  is  a  highly  non-linear 
function  of  the  impulsive  velocity  (V  ) .   The  critical  velocity 
is  very  dependent  on  the  value  of  y   used  as  shown  in  Figure  4. 

It  is  evident  from  Figure  4  that  the  displacement 
functions  at  T  =  T.  experience  larger  amplifications  than  do 
the  velocity  functions.   As  stated  by  Jones  and  Ahn  [13],  this 
suggests  that  shape  imperfections  might  exercise  a  much  more 
significant  influence  on  the  instability  of  spherical  shells 
than  variations  in  the  initial  velocity  field. 


* 
The  critical  harmonics  are  always  rounded  off  to  the  nearest 

integer. 
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V 
o 

(10  in/sec) 

Y 

-1 
sec 

Critical 
Mode 

A 
mn 

Critical 
Mode 
n 

B 
mn 

5 

40 

11 

1.21 

2 

0.12 

7.5 

40 

11 

1.47 

2 

0.18 

10 

40 

11 

1.85 

5 

0.24 

11 

40 

11 

2.05 

5 

0.26 

12 

40 

11 

2.28 

6 

0.29 

15 

40 

11 

3.17 

7 

0.40 

20 

40 

12 

5.81 

8 

0.70 

30 

40 

13 

21.6 

10 

2.38 

40 

40 

14 

88.6 

12 

9.07 

45 

40 

14 

332. 

13 

28.6 

5 

6500 

11 

1.60 

2 

0.19 

7.5 

6500 

11 

2.51 

6 

0.31 

10 

6500 

12 

4.22 

8 

0.50 

11 

6500 

12 

5.28 

8 

0.61 

12 

6500 

12 

6.64 

9 

0.75 

15 

6500 

12 

13.6 

10 

1.45 

20 

6500 

13 

48.8 

11 

4.81 

25 

6500 

14 

190. 

12 

17.4 

30 

6500 

15 

788. 

13 

68.4 

32.5 

6500 

15 

1644. 

14 

138. 

Table  3:   Critical  Mode  Numbers  and  Displacement  Function 

Amplitudes  Versus  Initial  Velocities  at  Differing 

Values  of  y   with  R  =  1.460  in.,  R/h  =  18.25, 

o0  =  46,000  psi  and  p  =  0.0002523  lb-sec2/in. u .  (Cont. ) 
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V 
o 

(10  in/sec) 

Y 

Critical 

A 
mn 

Critical 

-1 
sec 

Mode 
n 

Mode 
n 

B 
mn 

5 

1,288,000 

•  11 

1.99 

5 

0.25 

7.5 

1,288,000 

12 

3.85 

8 

0.48 

10 

1,288,000 

12 

8.26 

9 

0.96 

11 

1,288,000 

12 

11.4 

10 

1.28 

12 

1,288,000 

13 

15.9 

10 

1.76 

15 

1,288,000 

13 

45.7 

11 

4.73 

17.5 

1,288,000 

14 

116. 

12 

11.4 

20 

1,288,000 

14 

304. 

13 

28.6 

25 

1,288,000 

15 

2334. 

14 

205. 

5 

900,000,000 

11 

2.24 

6 

0.29 

7.5 

900,000,000 

12 

4.89 

8 

0.61 

10 

900,000,000 

12 

11.7 

10 

1.35 

11 

900,000,000 

13 

17.5 

10 

1.98 

12 

900,000,000 

13 

25.7 

11 

2.81 

15 

900,000,000 

14 

89.4 

12 

9.19 

17.5 

900,000,000 

14 

274. 

13 

26.6 

20 

900,000,000 

15 

856. 

13 

80.3 

22.5 

900,000,000 

15 

2902. 

14 

252. 

Table  3:   Critical  Mode  Numbers  and  Displacement  Function 

Amplitudes  Versus  Initial  Velocities  at  Differing 

Values  of  y  with  R  =  1.460  in.,  R/h  =  18.25, 

a      =    46,000  psi  and    =  0.0002523  lb-sec2/in. 4 . 
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It  should  be  noted  that  the  present  theory  is  an 
approximate  one  since  the  linearized  from  of  the  constitutive 
equations  for  rate-sensitive  material  were  used.   In  addition 
to  the  various  simplifications  and  approximations  which  were 
introduced  in  the  theoretical  analysis,  the  influence  of 
transverse  shear  stresses  and  deformations  have  been 
disregarded.   Also,  the  influence  of  unloading  has  been 
neglected  so  that  the  effect  of  possible  plastic  behavior 
after  t  =  Tf  has  not  been  considered.   Moreover,  the  influence 
of  wave  disturbance  has  been  neglected.   However,  with  final 
times  (t.p)  being  in  the  neighborhood  of  10-30  usee,  which  can 
be  calculated  by  using  equation  (lie)  and  Table  2c,  then  one 
should  not  neglect  wave  phenomena.   As  an  example  materials 
1,  10,  20,  30,  40  have  values  of  tf  =  8.7,  13.2,  11.3,  29.1, 
and  15.9  ysec,  respectively. 
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Chapter  6 
Conclusions 

A  theoretical  investigation  has  been  undertaken  into  the 
dynamic  instability  of  complete  spherical  shells  which  are  made 
from  a  viscoplastic  material  and  loaded  impulsively. 
Unfortunately,  no  experimental  results  are  available  to  examine 
the  validity,  or  otherwise,  of  the  various  theoretical  predictions, 
However,  the  forms  of  the  expressions  for  the  critical  mode 
number  and  threshold  velocity  are  similar  to  those  developed 
for  cylindrical  shells  and  spherical  shells,  except  that  the 
magnitude  of  the  numerical  coefficients  are  different.   The 
threshold  velocities  for  viscoplastic  spherical  shells  came  out 
higher  than  that  for  rigid  plastic  spherical  shells  with  the 
same  R/h  ratios  and  material  parameters.   The  viscoplastic  theory 
predicts  critical  mode  numbers  for  spherical  shells  which  are 
similar  to  those  observed  experimentally  on  cylindrical  shells 
with  10  <  R/h  <  30  but  the  correlation  was  better  at  the  lower 
range  the  critical  mode  numbers  were  only  one  to  three  less 
than  the  rigid-plastic  prediction  where  as  the  higher  range 
produced  mode  numbers  of  four  to  eight  less  than  the  rigid 
plastic  prediction. 
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Chapter  7 
Recommendations 

To  further  the  investigation  of  viscoplastic  effects,  the 
next  logical  step  would  be  to  look  into  the  non-linear 
constitutive  equations.   This  may  result  in  a  sufficient  enough 
change  so  as  to  correct  some  of  the  problems  of  the  threshold 
velocities. 
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Appendix  A 
Constitutive  Equations  of  Viscoplastic  Spherical  Shells 


The  constitutive  equations  for  cylindrical  shells  which 
are  made  from  a  rigid  viscoplastic  material  are  developed  in 
references  [8],  [10],  [11],  and  [18].   This  appendix  contains 
a  brief  derivation  of  constitutive  equations  for  viscoplastic 
spherical  shells. 

Consider  a  special  linearized  case  of  the  constitutive 
equations  for  a  rate-sensitive  plastic  material  as  derived  by 
Perzyna  [18] 


dF 


ID 


=  2ky  F 


(1.1) 


3a. 


ID 


where  e.   .  and  a. .  are  components  of  the  plastic  strain  rate 

and  stress  tensors,  k  is  the  yield  stress  in  simple  shear  and 

Y  is  a  physical  constant  of  the  material  later  termed  the 

viscosity  constant.  y   is  obtained  by  uniaxial  tests  as  described 

later.   In  the  above  equation  the  function  F,  a  function  of 

strain  rate  and  stress,  is  defined  by  F  =  /jT/k  -  1  where 

J0  =  1/2  S. .  S. .  is  the  second  invarient  of  the  stress  deviator 
2        13  xj 


ow 


with  components  S . . .   In  accordance  with  the  Levy-Mises  f 1 

rule  for  rigid  perfectly  plastic  material,  as  y   -*■   °°,  J-  =  k2 

where  k  =  a  //3. 
o 

The  results  of  uniaxial  stress  experiments  produced  the 
Cowper  and  Symonds  [19]  power-type  relation 


o      =    o 
e     o 


1  + 


e 


I  Y  J 


(1.2) 
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which  give  the  non-linearized  form  of  function  F. 
(1.2)  we  get 


Inverting 


ee  =  Y 


a 

-£  -  1 


(1.3a) 


3  h 

where  o      and  £   are  the  equivalent  stress  [—  S . .  S. .]   and 
e       e  ^  2   ID   ij 

2  •    Jj 

strain  rate  [—  £  .  .  e. .]   respectively  and  o      is  the  static 

3  13   ij  o 

uniaxial  yield  stress  of  the  material  and  y  an<3-  5  are  constants 
obtained  experimentally  for  a  given  material 


2  2        2  % 

a      =  [  cr ., ,   -*-a,..  a00+a00   +  3  a .  ~  I 
e      11      11   22     22       12 


2   r-   2 
e    /3 .11 


2    ,.  •   2  ^ 
11  £22  +  £22   +  ,£12  * 


(1.3b) 
(1.3c) 


Taking  the  linearized  form  of  equation  (1.2)  we  get 


0=0 
e    o 


1  + 


Y 


(1.4) 


where  y  is  redefine  as  explained  in  Appendix  E. 
In  the  case  of  plane  stress  (a-,-,  =  a„.,  =  a 
equation  (1.1)  yields  the  following  equations. 


33 


=  0) 


11 


Y 

3 


VT2 


-     1 


2a 


11 


22 


(1.5a) 


22 
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/j-2 


-    1 


2a22  -  axl 


(1.5b) 


E12  +  E21  =  2Y 


-  1 


12 


(1.5c) 
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where 


2    3   e 


k  =  a  //3. 


Inverting  equations  (1.5a-c)  produces, 


(1.5d-e) 


'11 


22 


12 


-r  (2eu  +  e22) 


(2e22    +  £]L1) 


E      +   F 

12     21 


(1.6a) 
(1.6b) 
(1.6c) 


where 


Y*  =  yd  -  k//5T)  =  yd  -  o/o). 


(1.6d) 


2'     '  '      o'  "  e' 
Now,  assume  that  the  strain  rates  can  be  represented  as  follows: 


ell  =  el  +  2  Kl 


e22    e2  +  Z  K2 


E12  +  E21  "  S12  +  Z  K12 


(1.7a) 
d.7b) 
(1.7c) 


Considering  that  y*  is  not  constant  and  in  fact  is  a  function  of 
z,  one  assumes  that  the  variation  of  y*  with  z  is  second  order, 

a  perturbation  quantity.   As  shown  in  appendix  B  y*  is  dependent 

•  •     jt_     .  i        .i 

on  e   and  as  shown  in  Figure  5,  e,  =  en  +  c.  where  e,  is  a 
e  ^111        1 

perturbation  quantity  dependent  on  z,.   Therefore,  y*  can  be 

taken  as  a  function  of  e   and  not  a  function  of  z.   Now  the 

e 

constitutive  equations  can  be  integrated  out  as 

h/2 

kh 


N, 


a  dz 


(2e1  +  e2) 


(1.8a) 


-h/2 
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h/2 


N2  = 


kh 


°2dz         y*     —2  ■  -1' 


(2e„  +  en) 


-h/2 
h/2 


N 


12 


rt   ,    kh  eI2 

a, -dz  =  —  

12      Y    2 


-h/2 


(1.8b) 


(1.8c) 


h/2 


Ml  = 


an  zdz  = (2k..  +  k0) 


12Y 


*  »-"l    ~2 


-h/2 


(1.8d) 


h/2 


M   = 
M2 


a~zdz 


kh 


12Y 


-  (2<2  +  ie  ) 


-h/2 


(1.8e) 


h/2 


M      = 

™12 


a,  2Z(^Z 


kh3   K12 


12Y*    2 


-h/2 


(1.8f) 
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R+h/2 


R-h/2 


Figure  5:   Variation  of  Strain  Rate  Across 
Thickness  of  Section. 
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Appendix  B 
Variations  of  y   and  a 


Equation  (1.4)  from  Appendix  A  is  shown  in  figure  6  along 
with  a  straight  line  tangent  to  the  Cowper  and  Symonds 
power  type  relation  as  suggested  by  Florence  [7] .   The  equation 
Florence  suggested  is 


a      =   oni    +   ye 
e     01      e 


(2.1) 


where  y  is  the  slope  of  the  tangent  line,  tangent  at  the 

•  * 

maximum  initial  strain  rate  e  ,  and  cA,  is  the  intercept  of 

e       01  c 

this  tangent  line  at  static  test  strain  rate.   The  slope  of 

the  tangent  line  is  obtained  by  differentiations  of  equation  (1.2) 

and  evaluating  the  results  at  e 

e 


da 


de 


y  = 


e  =  e 
e    e 


a- 

1      e 


(2.2) 


The  value  of  a  ,  is  now  obtained  by  equation  equation  (2.1)  and 


•  * 


(1.2)  evaluated  at  e 


01    o 


1  + 

e 

\ 

r 
i  - 

i 

[y  j 

. 

6 

J 



_ 

(2.3) 


Putting  equation  (2.1)  into  the  form  of  equation  (1.4) 


e     01 


1  + 


a01/y 


(2.4) 
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where  now  a.,  and  o~^/\i   compare  with  a   and  y    from  equation  (1.4) 
respectively.   Letting  a0n/v  =  Yi 


6.  -K 

•  l     e 


Y 


1/6 


1  - 


(2.5) 


and 


ae  =  a01 


1  +  -£ 
Yi 


(2.6) 


A  replacement  of  a   and  y   with  a_,  and  y,  respectively  in  the 
equations  would  give  the  identical  solution  with  the  result 
being  the  use  of  the  tangent  line  approximation. 

The  values  of  the  material  constants  are  6=5  and 
y  =  40.4  sec  '  for  hot  rolled  mild  steel  [20]  and  6=4  and 
Y  =  6500  sec    for  6061  T6  aluminum  [21] .   For  thin  spherical 
shells  of  radius  R  the  initial  strain  rate  is  computed  as 

.  £ 

e   =  2V  /R,  where  V   is  the  uniform  initial  velocity. 
e     o  o 


From  Chapter  2  equation  (5b) 

Y*  =  Yd  -  VCe} 

where  k   =   a   //3      and      /J0    =   o//3. 

o  2     e 


(2.7a) 
(2.7b-c) 


Solving  equation  (1.4)  for  cr  /ae  gives 


e   +  Y 

e    ' 


(2.8) 


Since  1/y   is  used  more  predominantly,  substituting  (2.8)  into 
(2.7a)  and  inverting  produces 


*      • 

Y     Y    ee 


(2.9) 
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(k  /c  2) 

g   ra 


1/a 

a    =  a 
e        o 

1    + 

.  Y  . 

e 


ee  (sec   ) 


Figure  6:   Cowper  and  Symonds  Pov/er  Law  and  Linear 
Approximation  to  it  [11] . 
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Appendix  C 
Computer  Proaram  Usina  Florence  Method 


The  following  program  was  used  to  solve  for  A  n(£)  and 
B   (£)/  equations  (36a)  and  (36b),  by  use  of  Mathematics 
Library  Subroutine  called  HYPER  which  is  a  subroutine  to  solve 
Kummer ' s  function  of  the  first  kind.   Problems  were  encountered 
when  taking  the  difference  between  functions  as  this  was 
eliminating  the  first  five  or  six  significant  figures. 
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Appendix  D 
Computer  Program  Using  Numerical  Method 

The  following  program  solved  the  basic  equation  (34a) 

for  A   (£)  and  B   (£)  by  use  of  Mathematics  Library  Subroutine 
mn  mn      J  J 

called  DQSYK  which  is  a  subroutine  to  solve  up  to  fourth  order 
simultaneous  differential  equations  with  initial  conditions. 
This  program  used  equation  (34d) . 
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Appendix  E 
Critical  Impulse 


It  is  straightforward  to  show  when  A   >>  1  that  X      can 

n  n 

be  eliminated  from  equations  (34c)  and  (34e)  and  using 
equation  (24b)  to  give 


pV  2    8     R  2  h/R 

o  n  ' 


/3  (a  /a  )  2 
v  o'    e 


24/3  Q(2) 
n 
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1/2 


(5.1) 


where  o   /a   is  defined  in  equation  (30).   As  shown  in  Table  4 
o   e 

the  values  of  a  /a   lie  within  the  range  0.98  <  a    >'o      <  1.16 

o   e  o   e 

(2) 
with  a  medium  value  o   /a      -    1.07.   Also  the  values  of  O     lie 

o   e  ~n 

(2) 
within  the  range  1.27  <  0     <  4.42  with  a  medium  value  of  2.84 

It  is  straightforward  to  show  when  using  equation  (34c)  with 

n  >>  1,  that  the  critical  harmonic 

h 

(2) 

is  associated  with  these  values  of  a  /a   and  Q    . 

o   e      n 


n  =  3.173 (R/h) 


(5.2) 


The  threshold  of  dynamic  plastic  instability  could  be  said 
to  occur  when  A  n(0)  =  X,  where  X  is  an  arbitrary  positive 
number.   If  X  =  100,  then  the  initial  displacement  imperfections 
in  the  critical  mode  are  magnified  by  a  factor  of  one  hundred. 
This  amplification  is  achieved  when  R   -1.7  according  to 
Florence  [3]  asymptotic  formula  for  equation  (36a) 
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Amn(0)  =  exp 


-(A  +  n  ) 

n    n 


(A  +l+b  )M(S2  )+fi  M'  (fi  ) 
n    n    n    n    n 


1  +  b 


n 


(5.3) 


where  the  right  hand  quantities  are  defined  in  equations 
(37a-h) . 


,  V° 
(103in/sec) 

(sec   ) 

a  /a 
o7  e 

Critical 
Harmonic 

Q(2) 
^n 

5 

40 

1.00 

11 

1.27 

30 

40 

0.98 

13 

2.23 

5 

6500 

1.06 

11 

1.35 

30 

6500 

1.03 

15 

3.99 

5 

1,288,000 

1.11 

11 

1.41 

25 

1,288,000 

1.10 

15 

4.26 

5 

900,000,000 

1.16 

11 

1.47 

22.5 

900,000,000 

1.14 

15 

4.42 

Table  4: 


Variation  of  a   /o      and  O 
o'     e      n 


(2) 


to  various  values  of 


and  y  with  R  =  1.460  in.,  R/h  =  18.25  a      =  46,000  psi 
and  p  =  0.0002523  lb-sec2/in. h  . 


Finally  equation  (5.1)  predicts  a  non-dimensionalized  critical 
or  threshold  velocity 


PV 


1.11  h/R 


(5.4) 


This  asymptotic  formula  gave  agreement  with  numerical  values 
at  £  =  0  up  to  4  significant  figures. 
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This  value  does  not  compare  favorably  with  Jones1  [13]  in  which 
his  coefficient  is  17.8,  since  the  velocities  required  should 
be  higher  than  in  the  in  rigid-plastic  case. 
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